function x = LxbSol(L,b,n)
%UnitLxbSol 解单位下三角阵方程
%输出结果x储存在b中
for j = 1:n-1
    b(j) = b(j)/L(j,j);
    b(j+1:n) = b(j+1:n) - L(j+1:n,j)*b(j); 
end
b(n) = b(n)/L(n,n);
x = b;
end

